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FOREWORD 


Under contract to Union Carbide | the Jet Propulsion Laboratory 
has undertaken the design and construction of a temperature controlled 
device for the freezing and thawing of red blood cells. Stringent 
temperature control is obtained from electrical circuitry used on the 
Viking Mars Mission for thermal control of spacecraft instruments. 

This control methodology allows precision freezing and thawing in a 
reproducible manner and at rates which may be optimal for maximum 
recovery of the red blood cells. 

This is the final report for the project. 

The development phase of the project is described in the Red 
Blood Cell Freezer Task Interim Report, JPL internal document number 
5030-134, dated October 1, 1977. 
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ABSTRACT 


A system has been constructed which allows programmable 
temperature-time control for a 5cc sample volume of an arbitrary 
biological material. It also measures the parameters necessary for 
the determination of the sample volume specific heat and thermal 
conductivity as a function of temperature^ and provides a detailed 
measurement of the temperature during phase change and a means of 
calculating the heat of the phase change. 

Steady-state and dynamic temperature control is obtained by 
supplying heat to the sample volume through resistive elements 
constructed as an integral part of the sample container. For cooling 
purposes, this container is totally immersed into a cold heat sink. 
Sample volume thermodynamic property data are obtained by measurements 
of heater power and heat flux through the container walls. Using a 
mixture of dry ice and alcohol at -79®C, the sample volume can be 
controlled from +40®C to -60^C at rates from steady state to 
± 65^C/min. Steady-state temperature precision is better than 
0.2^C while the dynamic capability depends on the temperature rate 
of change as well as the thermal mass of both the sample and the 
container. 

Preliminary tests, conducted with distilled water, saline, and 
red blood cells, showed the feasibility for using the device for 
extremely accurate temperature control. The performance capabilities 
of the system allow the application of an increased dimension of 
thermal control to the study of freeze storage of biological 
materials. The system’s ability to determine the temperature 
dependent thermodynamic properties allows a greater understanding of 
the thermodynamics of freeze damage to biological cells. The unique 
capability in determining phase change and measuring its heat, may 
lead to a new consideration of the roles of the electrolytes through a 
more thorough understanding of the solution eutectics and their 
resulting osmotic pressure histories. 
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SECTION I 

INTRODUCTION AND BACKGROUND 


Numerous attempts have been made to freese biological materials 
and show that viability is maintained upon thawing. In some cases, 
satisfactory freeze preservation has been possible when cryophylactic 
agents, such as dimethyl sulfoxide (DMSO) or glycerol, have been added 
to materials before freezing. The exact manner by which these 
materials afford protection is not known, but by understanding the 

various temperature related factors which cause cryogenic destruction, | 

some understanding of the effects of freeze/thaw protocols may be | 

possible. ! 

i 

The problems arising during the freezing process include: (1) > 

intracellular and extracellular icc formation, (2) high solute •’ 

concentration, (3) pH changes, (4) dehydration, ard (5) others. It is 
known however, that the freeze/ thaw rates and storage temperature 
affect both the character of ice formation and other variables which 
indirectly influence the survivability of frozen materials. 

Temperature rates which are either too slow or too fast are known to 
cause increased cellular destruction. 

The use of frozen thawed human erythrocytes for transfusion is 
becoming increasingly important for many reasons, including, for one, 
the convenience of having rare blood types available when needed. The 
fact that frozen blood is not even more widespread in its use can be 
attributed to its relatively high cost in relation to the unfrozen 
product. The high cost is largely because large quantities of 
glycerol must be added as a cryoprotectant and then later removed. 

The removal of glycerol is a complex, time consuming, and costly 
procedure . 

This work describes a device which allows the study of a new 
approach to the freeze preservation of human erythrocytes and which, 
if successful, will have a significant impact on the frozen blood 
industry. The approach utilizes technology developed at the Jet 
Propulsion Laboratory (JPL) to allow precise control of instrument 
temperature aboard spacecraft. With this precise control, it may be 
possible to understand the role of temperature for minimizing freeze 
damage, thereby minimizing the costs for low-temperature preserva- 
tion. Application of the JPL-designed temperature control system to 
the problem of red cell freezing may provide solutions to these and 
other problems of freeze preservation. 
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SECTION II 


TECHNICAL DISCUSSION 


A. OVERALL SYSTEM DESCRIPTION 

1 . Concept 

In the initial phase of this work, a basic design concept was 
adopted which allowed the red blood cell container to also serve as 
the primary means for heating the cell mass while immersed within a 
cold environment. Such an approach maximized utilization of the 
temperature control available from the bridge circuit proposed for use 
and discussed in Appendix A of this report. This design concept is 
summarized in Figure 1. The desired freezing and/or thawing protocol 
is established in the computer. The computer measures the heat 
dissipation in the heater and adjusts it by sensing the temperature of 
the cell mass or other control point and referencing this against the 
desired temperature for a particular protocol. The heat of phase 
change as well as the thermal conductivity and specific heat of the 
cell mass arc also determined from temperatures and heat dissipated in 
the container. These parameters are necessary for thermal dimensional 
scaling to larger volumes of cells which may be of commercial interest 
in the future as well as for the normalization of the thermodynamic 
behavior of the red blood cell samples obtained from various donors. 


2. Approach 

The final system block diagram is given as Figure 2. The 
function of the various blocks of this configuration are discussed in 
this section. 


a. Temperature Controller . The controller sets tne 
temperature levels for the system. It is discussed in Section II-C-1. 


b. IDaC . The IDAC provides automatic operation of the 

controller, transmits digital data to the tape deck, and produces a 
video display in engineering units for any selected input. The IDAC 
is described in Section II-D of this report. 


c. Freeze/Thaw container - Model III , The configuration of 
the container used for the final experiments in this report is 
described in detail in Section II-B-2 of this report. Samples of 
water, saline solutions, red cells, or other liquids can be placed in 
the container for experimental purposes. The sample volume is 5 cc. 

A thermocouple is located in the center of the sample. The container 
walls are hard anodized aluminum separated by a Vespel rectangular 
annulus. The walls contain the temperature sensors and heaters which 
provide the total thermodynamic information for the system. 


<!• Digital Voltacttr . Any digital voltmeter (OVM) with an 
input impedance greater than 100 kli and a resolution of least 1 /iV 
can be used with this system. The DVN is used to monitor the various 
temperature sensors during operation or calibration, and to monitor 
the sensor excitations and controller sensitivity during system set up. 


e. Heater Power Supply . Any dc power supply that can provide 
150 Vdc at 10 A can be used with this system. Typical operation is 
100 Vac and 7.5 A at full power ir the alcohol-frosen CO 2 heatsink. 

The dc supply provides dc power to the temperature controller which in 
turn provides the pulse width modulated power to the heaters in the 
container walls. 


X-Y plotters . The 2 pen plotters are used to graph any 
combination of sensor and heater signals desired. The plotters 
provide real-time demonstrations of the performance of the system. 


g. Tape Deck and Video display . The tape deck stores the 
complete thermal history of each test. The video display is used for 
monitoring up to four channels of real-time test data. These are 
discussed in detail in the IDAC section of the report. 


B. FREEZER DESCRIPTION 

1. Early Development 

Two concepts wer>:! considered in the initial design phase. The 
first, based on a rotating cylindrical tube container, is discussed in 
Appendix B. This approach was determined to be inappropriate for the 
basic studies initially required. In the second concept, a container 
having parallel flat plates separated by a guard material would be 
used. The entire heat transfer between the cell mass and the heat 
sink would occur across the flat plates. These plates would also 
contain the necessary heaters and sensors. This approach was chosen 
for the study because of its simplicity and ease of analysis. 

A schematic of an early model of this configuration, Model I 
shown in Figure 3. was built for the primary purpose of exercising the 
temperature control circuit. Model I consisted of a pair of plane 
parallel 3 x 6-in. walls held apart by an 0.08 in. aluminum spacer and 
enclosed in a plexiglass rim. Each w^xll contained a pair of O.OAO-in. 
beryllium copper plates which sandwich(‘d an 0.080 in. Electrofilm 
patch heater. The heater and plates were bonded together by RTV 60 
cement. 


The temperatures for Model I were monitored by a pair of 0.5-in. 
square nickel sensors placed at the center of the wall between the 
inner plate and the patch heater. The interior volume of 31 cc was 
used to contain distilled water for the initial experiments. 
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The Actual thermal perfoimance of Model although inadequate for 
blood freezing aCudieat provided data which suggested design nodif icatioi.a 
for the control electronics and the subsequent contaioei configurations* 

These changes led to the design of Model II ^ shown in Figure 4« The 
primary purpose of this model was to test the thermal control capability of 
specially designed heater and sensor elements* 

In Model II » the cell mass is contained between a pair of 4 x 6 in* 
plane parallel wallSi separated and insulated from each other by a styrofoam 
rim and an 0.080*in* Vespel SP21 Spacer which served as a guard volume for 
the cell mass* Vespel has thermodynamic properties similar to packed red 
cells. It introduced no noticeable phase change during the freeze/thaw 
protocols being considered* The wall itself is composed of four elements, 
bonded to each other with Micro-measurements M-Bond 610* This bond was 
successfully tested under extreme thermal stresses by immersion of various 
combinations of bonded surfttces into liquid nitrogen from room temperature* 

No separations were experienced even though the bonded materials had 
considerably different thermal coefficients of expansion* 

The interior element of the wall was a 0.032 in* anodized aluminum 
plate whose primary function was to provide a high thermal conductivity path 
be^wee^ the cell mass and the control sensor, the adjacent element of the 
wall. The control sensor was a 0.0005 in. nickel grid deposited on a 0*001 
in. Kapton sheet* This sensor provided an integrated measure of the cell 
mass surface temperature on each interior wall of the container* The next 
element of the wall is physically the same* 

The outer element of the wall was a 0.090 in* anodized aluminum plate 
whose primary function was strength with a high thermal conductance path to 
the heat sink. Because of the reduced thermal mass of Model II over Model I, 
Model II was used to indicate the true order of magnitude of the temperature 
control available from the pulse width modulated heater control system* 
Experiments with Model II were conducted to optimize the heat transfer 
performance of the container* 

Model IIIA, shown in Figure 5, was also considered in concept but never 
constructed. It had the potential of exhibiting a maximum thermal coupling 
between the cell mass, its surface sensor and the heater* This would have 
been accomplished by depositing both the nichrome heater and the nickel 
sensor on a single sheet of 0*001 in. Kapton* Electrical insulation between 
the heater and sensor and between the sensor and the cell mass would be 
obtained by a thin deposited layer of silicon monoxide. The remainder of the 
wall elements would be the same as described for Model III* 

Although Model IIIA would have required some additional development 
work, it would have had other important advantages* Its basic construction 
would simulate the presence of a plastic bag by the inner Kapton sheet. This 
simulation could be important in a developmental sei.se in that the. use of a 
plastic blood bag for freezing and storage is commercially attraci.xve* 
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Final Freezer Configuration 


a. Model III Description . Model III was designed and 
constructed as a container for red blood cells having sensors to 
accurately measure temperatures, temperature gradients, heat flow and 
state changes during the freeze/thaw process. It is operable in the 
temperature range from room temperature to -60^C when used with a 
heat sink of -79°C. The rate of cooling/heating can be controlled 
either from the data sensed by a thermocouple in the cell mass or from 
sensors in contact with the cell mass surface. The latter method was 
chosen for these experiments. The directional heat flow is measured 
by sensing elements within each wall. This heat flow measurement 
provides data to determine the thermodynamic properties of the cell 
mass and to measure the heats of phase change. 


b. Description and Explanation . Figure 6 is a schematic 
representation of the Model III container. The Model III contained a 
total of 5 ml of cell volume with 112 expansion allowable for volume 
change during freezing. The Model III design was developed from the 
results of earlier tests. 

Each wall of the container is composed of 9 individual elements 
as shown in Figure 6 . Each element is attached to the adjacent 
element with Micromeasurements AE-15 bond. The two walls sandwich the 
Vespel rectangular annulus (V) allowing a 2mm spacing of the walls and 
the formation of the 6.5cc sample volume. Material is transferred to 
and from the sample volume by syringe which can be inserted into ports 
drilled through the Vespel frame. A hydrophobic grease is placed 
between the Vespel and the walls to eliminate leaks and a spring 
provides the necessary tension to hold both walls and Vespel in place. 

All the plates shown in each wall are hard anodized aluminum. 

The anodized coating provides electrical insulation between the plates 
and the sensors or heaters, and provides an acceptable surface for 
contact with the sample mass. The positioning of the aluminum plates 
between the heaters and sensors provides the necessary electrostatic 
shielding to prevent the pulsed heater power from coupling into the 
sensors. 

The dimensions of the wall elements are 2.625 x 3.35 in. The 
sensors shown in Figure 7 are 0.0005 in. nickel deposited on 0.001 in. 
Kapton. Both heaters and sensors are bifilar wound to reduce 
electrostatic coupling. The total heater resistance is 25 ohms and 
the total sensor resistance is 94 ohms. The internal aluminum plates 
are all 0.020 in. thick and the two outer plates, R 5 ^ and R 5 g, are 
0.080 in. thick. 

Sensors and measure the surface temperature of the 
sample mass. The temperature differential across a portion of the 
wall is made from the S 2 a ~ 83 ^ and the S 2 g " 835 sensors. 

These later temperature differentials are important in calculating the 
heat flux through each wall. This heat flux is required for 
determination of the thermodynamic properties of the sample mass. 
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Fiuure 8 shows anotlu»r viow of tho coatsiuor walls and 
associatos oach sonsor and luviLor with its appropriate temperature or 
heater symbol used in this text* For example, the temperature 
measured by sensor power generated by heater 

Ha is Pa. 

Fip,ure 9 is a photograph illustrating a side view of the 
parallel plate oonstrvic tion ol one half of the container and showing 
the Vos pel frame on top* It also shows the heater and sensor leads 
coming from the wall. 

Figure 10 is a photograph of one half of the sample container. 

The Vespel frame is seen as the d^irk rectangular annulus. The heater 
and sensor element wiring from the wall uses plastic covered shielded 
twisted pair conductors. 

Figure 11 shows the outer platt' of one si vie, together with the 
styrofv^iun edge insulators mounted in the lower half of the 
i ns t rumen t a t i on f r ame . 

Figure 12 shows the i nst rumentat i v>n fixture and illustrates the 
oxpanvlv'x spring mattM*ial whic!\ exerts iiiv'reasing prv'*ssure on the two 
halves of the container stav'k as t\w temperature is loweianl. This 
photograph alsv) shv>ws t lu' screw mv'ans fvu* jviini ng the two halves of 
the instrumentation fixture and illustrates the coolant flow apertures. 

c. Tvj>ical Mv^del^ 1 1 1^ Tt'S t batji. Figure 13 shows a typical 
temperature-time history for t ret'i'.ing anvl thawing of disti lled water. 
lVv> cov>ldv^wn xates weri' usevi: Hat v' I, 10 ^V/min, and Rate 2, 5 
^C/min. Curvt's AA ativl 88 ate ihc cot i't^sponding power curves for 
lu*ater As the sample^ UlistiMcd watv'r) is cvu>letl, it initiates 

a phase chaugt' at a temperatuiv* Inflow its uinaual tree/.ing pv>iut 
( supt'rcool i ng ^ ♦ The v'xt>t lu'uni c trt'er.ing pi ocess causes the 
temperature to iucieast' la Mu' diroctivm of the theoretical tree/Mng 
pv>int, O^V*. AfttM' til*' exotlumnic triu'iu' prv^cess has hocu completed, 
continued heat removal t rom the sample causes the tempo r iture to lirop 
rapivUv tv^ tlie ramp 1 ovuM M'ud ot Phase tUiange^* Further coolings 
continues at th<‘ ramp tate to tMthei' the next phase chiiuge point or to 
the terminat i t>n v>l the cooling iMiup. 

It should he nott'd Mi at the curvu's A and K are sinisoil by the 
t he ruK'tv' on p 1 v' . it t hv* vUMittM* of t hi' sample m iSvS. Had the control 

sensor S| heen plv>ftod, it woiiKl show a straight line through the 
phase change indicating pia'ci so temperatun' coi\t rol ot the actual 
coittiiM point. This characteristic is demons trated in the test data 
to be discussed. 


c. 


CONTROL SYSTEM DESCRIPTION 


1. Temperature Controller 

The data system temperature control block diagram is shovm in 
Figure 14. The temperature controller provides sensor excitation and 
pulsed heater power to the container, and digital information of 
controller set point temperature reference and the analog signals of 
the various nickel temperature sensors to the Integrated Data 
Acquisition and Control System (IDAC). Other analog outputs are 
provided for the X-Y plotters. Any of the freezer temperatures and 
heater powers can be plotted in real time. The controller provides 
the freeze/thaw control ramps and the rates at which these ramps are 
accomplished. The controller starts, stops and holds at manually 
selected control points or the programmed control points provided by 
the IDAC. All functions of the controller can be set manually or can 
be programmed in the IDAC. The IDAC then conducts a complete 
freeze/thaw experiment. The temperature controller circuitry is 
described in the following Sections. 

Figures 15 and 16 are front views of the controller chassis with 
and without the front chassis plate. Figure 17 is the back view of 
the chassis. A schematic showing the interconnection of the 
controller card circuit is shown in Figure 18. These card circuits 
will be discussed in the following paragraph. The temperature 
controller provides for both manual or programmed operation. 

a. Programmed Operation . For programmed operation the front 
panel run switch is in the hold position. Initially, it is desirable 
to have the freeze/thaw container at room temperature. The set point 
counter is then adjusted by the switches on the front panel until a 
proper setting is indicated by the output lights. The enable switch 
that provides a dc voltage to the power drivers, can then be turned on 
with appropriate change in the sample container temperature. The 
sample container is then placed in the coolant. The heater power 
immediately increases to the proper level to maintain the set 
tomperature. The IDAC program then sets the proper temperature rate 
and start temperature. The experiment program begins and the system 
accomplishes all of the programmed temperature rates, starts, stops, 
and holds. 

The control system utilizes sets of composite subprograms as 
building blocks to develop a total protocol. Each composite 
subprogram consists of (1) start temperature, (2) hold time at the 
start temperature, (3) temperature ramp rate, (4) final temperature, 
and (5) hold time at the final temperature. An example of this 
composite subprogram is shown in Figure 18a. 

The control system is capable of storing and managing up to 15 
of these composite subprograms per test protocol. Should more than 
this be required, the system can be placed into a manual hold mode 
while a new protocol is set in the computer. 


Should the start temperature for a given test not be equal to 
the current initial temperature of the system, the controller will 
take the system to that start temperature under the maximum rate 
possible. After holding the system at its final temperature for the 
appropriate time, the controller will automatically retain this same 
temperature for an additional indefinite time period or until it 
receives a manual or programmed command to do otherwise. 

b. Manual Operation . For manual operation the run switch is 
in the hold position and the set point counter is adjusted in the same 
manner as for the programmed operation. The sample container is 
placed in the coolant. The rate counter is set to the desired rate. 
The enable switch is positioned to on. To begin a cooling ramp, the 
run switch is placed in the down position. To stop and hold, the run 
switch is placed in hold position. The rate can be changed as 
desired. For thawing, the run switch is placed in the up position. 
Care should be exercised in manual thawing to prevent overheating the 
sample container. For manual operation, the sample temperature must 
be monitored continuously to achieve the desired start and stop points 

c. Front Panel Controls . The operating controls (Fig. 19) 
are front panel push button switches and toggle switches, or external 
voltage signals from the data acquisition system, that can set the 
ramp rate divider register and control the initial and final values of 
a temperature-time history ramp. The data acquisition system monitors 
the temperature bridge outputs and has an operating control program 
that provides the control logic for sequencing a series of ramps and 
holds . 


d. Up/Down Counter and Ramp Rate Clock . An up/down counter 
provides the 16-bit binary signals to control the D/A converter (see 
Figs. 19, 20 and 25). It counts up, down, or stop under control of 
external voltage or switch contact inputs. The count rate is 
controlled by a clock. The ramp clock rate is set by another digital 
binary counter that divides the source clock by numbers from 1 to 63 
in unit steps. The division ratio is controlled by an externally set 
register. This register is monitored and set by the external data 
acquisition system prior to the start of any voltage ramp out of the 
D/A converter, 

G. Comparator . This circuit (Fig, 21) is a high-gain, 
low-drift operational amplifier, LM108A, operating as a voltage 
comparator. One input is connected to the filtered output of the 
temperature bridge; the other input is connected to the output voltage 
of the controlled reference digital to analog converter plus the 
proportional band triangle wave signal (Fig. 24), The comparator 
output is on or off and provides drive to the switch mode power driver 

Power Driver . The power driver (Figs. 18 and 21) is a 
switching mode power transistor circuit that operates on the ground 
end of the heaters and completes the circuit through the heaters from 
the dc power supply. On-off switching times are in the microsecond 
range. The switching frequency is generated by the proportional band 
oscillator (400 Hz), and the duty cycle is controlled by the 
comparator output over the range of 0 to 100%. The power driver 
circuits also include 
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attenuators that provide conditioning of duty cycle signals and heater 
dc power for monitor by the exernal data acquisition system. The dc 
power for the heaters is provided by an adjustable, regulated dc power 
supply. The power supply must be capable of 150 Vdc at 10 A. 

g. Bridge Supply . The bridge voltage supply (Fig. 22 and 
Fig. 23) consists of a buffered set of stable inverting operational 
amplifiers that provides plus and minus 6.4 V dc for the temperature 
sensors. The reference source for the supply is the voltage reference 
output of the D/A converter. Since this reference is also used as the 
reference for the variable voltage output of the D/A converter, the 
temperature bridge operates in a full bridge compensated mode. 

h. Digital to Analog Converter . The D/A converter (Fig. 22) 
is a stable monotonic, 16-bit binary, voltage output D/A module. 
(Analog Devices. DAC HR16B). It provides the voltage reference for 
the bridge supply, and the variable output voltage controlled by the 
16-bit binary signals from the up-down counter to the comparator 
reference. The comparator reference is generated by adding a chopper 
stabilized operational inverter (the D/A output voltage) and an 
adjustable triangle voltage (Fig. 24) for proportional band set. The 
comparator reference voltage ranges from 3,2 volts to 6.4 volts for a 
binary input ranging from OOOOOOg to 177777g, This allows a 
potential control range from +40^C to -80®C. 

i. Container Heaters . The heater elements are thin film 
nichrome, in a bifilar configuration on a Kapton substrate. The 
heaters are located in the walls on each side of the freezing 
container. The heater resistance is 25 ohms. A photograph of the 
heater is shown in Figure 7. 

j. Container Temperature Sensors . The temperature sensors 
are thin film nickel in a bifilar configuration on a Kapton 
substrate. The sensors are located in the walls on each side of the 
freezing container. The nominal resistance of the sensor, at 22°C, 
is 94 ohms. A photograph of the sensor is shown in Figure 7. 

The power supplies for the controller circuits are shown in 
Figure 26. 


D. INTEGRATED DATA ACQUISITION AND CONTROL (IDAC) SYSTEM 

DESCRIPTION 

The Integrated Data Acquisition and Control (iDAC) data system 
was designed in 1966 as a general purpose digital data system for 
JPL*s propulsion test programs. The IDAC includes real time dcta 
output and control capabilities. The analog input subsystem is 
typical of many digital systems in that it has 500 channels, 14-bit 
ADC, scan rates to 10 kHz, addressable gain, and channel selection. 
The system includes an 8K, 18-bit memory and 78 in/sec digital 
magnetic tape recorder. Real-time engineering data is provided by a 
printed paper tape, two video displays, and 16 digital 
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to analog converters. An engineering units capability includes 
ranging and scale factors for all standard transducers, themocouple 
linearization, and data output unit definition. The system is able to 
limit check 64 data channels in real time. The IDAC also has a 
sophisticated digital data input capability. The digital input 
capability was used to provide the automatic control of the freezer. 
The operator supplies the freeze/ thaw control parameters to the 
control system of the IDAC. The system temperature controller 
supplies a reference input to the digital system. 

This digital reference input is compared with the freezer wall 
temperatures and the freeze/thaw control parameters. The result of 
the comparison is transmitted to the temperature controller which in 
turn sets the heater power of the freezer. The IDAC sets the rates of 
freezing and thawing, the start and stop temperatures and the hold 
times at the predetermined set points. The signal flow of the IDAC is 
shown in Figure 27. 


E. SYSTEM CALIBRATION 

The nickel sensors were calibrated by the JPL Standards 
Laboratory for resistance vs temperature. The temperature range of 
the calibration was from 22^C to -bO^C and produced the following 
equation: 

R » 83.699105 + 0.4065255 T + 0.000373 t2 

Calibration in the freezing system was accomplished by submerging the 
container in a distilled water-ice bath. When the system reached 
equilibrium, the bridge excitation voltages were adjusted to produce a 
O^C output from the data system. The system was also checked by 
submerging the container in room temperature water. The T^ 
thermocouple was used as a reference standard for these measurements. 
This is a type E thermocouple and was calibrated by the JPL Standards 
Laboratory and found to be within *01Z of the laboratory standard from 
22<5C to -60OC. 

The heater resistances have a nominal value of 24.3 ohms at 
0®C. This value was adjusted by a computer program, using the 
nichrome temperature vs resistance table, to produce an accurate 
measurement of each heater power. This is necessary since the heater 
temperatures were essentially the same as those of the control 
sensors, and changed with the cooling and warming protocols. 


1. Data Printout 

The typical output data presented in the following figures are 
obtained through an acquisition process conducted by the IDAC and 
which results in the data being stored on magnetic tape. A special 
program, modified from existing software, was prepared to accept this 
data, provide the necessary calibration and produce three listings of 
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the data in appropriate units. Typical examples of these listings are 
shown in the following figures. The first column of each listing is 
the time during the test in hours, minutes, and seconds, (HH:MM:SS)« 

Figure 28 presents the basic temperature and power data taken at 
one second intervals during a portion of a test. The temperatures 
Ta 3 , Ta 2 i and T^i are on the outside wall of the container 
adjacent to the cooling bath, within the package wall and at the 
interior face adjacent to the sample mass respectively. The 
temperatures T 33 , Tg 2 > and Tg]^ are in similar positions on the 
opposite wall of the container. The temperature at the center of the 
sample mass is Tq. The temperature set as a voltage within the 
controller to define the required temperature for a given protocol is 
Trep* The power dissipated in the container wall heaters is given 
as ?A and Pg, 

Figure 29 presents data characteristic of the system 
performance. The terra (Tq-Tref) indicates the difference between 
the actual temperature as measured at the center of the sample mass 
and that set by the controller. This is a measure of the overall 
ability of the controller to influence the temperature of the least 
responsive area of the sample mass. The term indicates the effects of 
thermal inertia, and phase change, as well as sensor and thermocouple 
calibrations. The value of (Tq-Tai), is the gradient from the 
center of the sample mass to the surface of the sample mass on the (A) 
side. The term (T^-Tg]^) is similar except on the (B) side. These 
terms also reflect the sensor calibrations and are useful for 
indicating the velocity of the wave front during phase change as well 
as the thermal conductivity of the sample mass. 

The terms (Ta 2 ”Ta 3 ) and (Tg 2 ”Tg 3 ) are the differences in 
the temperature between a plane within the container wall and one at 
the exterior surface where the container wall comes into contact with 
the coolant. These terms provide a measure of the heat transfer 
through the container walls and establish the thermal conductance of 
the walls used to evaluate the thermodynamic properties of the sample 
mass. Approximate values of these thermal conductances are calculated 
and presented in the columns under the title Ka and Kg of Figure 
29. They are “effective** values which are true measurements of the 
exterior wall conductance only under conditions of steady state. 

Under transient conditions, their perturbations from the true value 
may reflect conditions of non-equilibrium due to thermal inertia or 
phase change. The final column, average power, is that calculated by 
taking one-half of the sum of the power from the (A) and (B) side 
heaters. This term is used for plotting the thermodynamic cycle 
presented by each protocol in the test data section of this report. 

Figure 30 presents data of general engineering interest. The 
terra (Ta 1 “Tref) indicates the control system* s ability to produce 
changes in an area of maximum sensitivity to the control heater, i.e,, 
at the sample mass/container wall interface. Typical data for steady 
state conditions are shown as curve (a) of Figure 32. 
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Because of the importance of making accurate measurements of the 
outer wall conductance, a second indication of the temperature 
gradient across the outer wall of the container is available by a 
differential measurement. These values are presented as "Measured^' 
(Ta 2 “Ta 3 ) and "Measured" (Tb 2 "T 33 >. They can be used for a 
more detailed analysis of the data. 

The term "Heater Source" (E3) is used in the calculation of the 
power dissipation in heaters (A) and (B) and is the square of the 
power source voltage. The adjacent column is the calculated 
resistance used for heaters (A) and (B) and which is obtained from the 
appropriate calibration curve as a function of temperature. The 
column labeled "Shorted Channel" is the effective temperature produced 
by an unused channel, shorted near the experimental freezer. It 
therefore is an indication of the noise from the cable, the 
multiplexer scanner, the multiplexer amplifier, and the A/D 
converter. These data are used to indicate the failure of any major 
part of the control and acquisition system which processes data in 
parallel. Noise levels indicated here can be directly added to the 
performance curves of Figure 31a but are negligible. 


2. Steady-State Operation 

Curve (a) of Figure 31a presents the difference between the 
steady-state thermocouple temperature as measured in the center of the 
sample mass (Tq), and that established by the controller (Tggp). 

The data was obtained for water but should be typical for any sample 
mass. Examination indicated that this error was consistent and could 
be removed from the measurements. An approximate correction can be 
made to curve (a) of Figure 31a with the following equation 

(Tq - Tj^EF^ corrected ~ “9.8271 x 10 ^ , 

+ 1.0654 X 10"2 Tset “ 1-344 x 10~1 
+ (To - Tref). 

With this correction, the maximum error is 0.2^C as shown in Figure 
31b. This correction is not applied to the test data presented in 
this report and consequently it will appear that a hold at a low 
temperature will be in error a few degrees. 

Curve (b) in Figure 31a shows the difference between the 
temperature measured at the B sensor (Tgj) and that set by the 
controller (TggF^* This is an indication of the control capability 
of the system and does not involve the sensor calibration. The curve 
is repeated on a magnified scale as curve (a) of Figure 32 which shows 
that there is less than O.IS^C error in steady state temperature 
control and this error is greatest at the higher temperatures. In the 
range below O^C, a more important area for thermal control, this 
error is less than O.IO^C. 

Curve (b) of Figure 32 shows the difference between the temper- 
ature requested of the control system (TgET^ liiat interpreted by 

the control system as the desired set point (Tj^gF^* difference 
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is an indication of the inability of the control system to read and 
establish the proper set point* Since the error has an approximate dc 
level around ♦0.05®C, it may be possible to remove some of that 
error through a bias of the temperature submitted to the controller* 

If the error of curve (b) is due only to the dead band of the 
controller) then adjustment may not be possible* 

Curve (c) of Figure 31a is the net error between the desired 
steady-state temperature that measured at surfaces of the 

sample mass (T^i) and (Tg^)* The (A) face temperatures are indi- 
cated by the open triangles and the (B) face by the solid ones* Curve 
(c) is the maximum expected error of the control system in producing a 
repeatable steady state temperature* If a bias of the (Igg-p) can be 
made to correct for the error sho^^ in Curve (b), then the error 
indicated in Curve (c) may also be reduced by an equivalent amount* 

In addition, the characteristic syrtem ao: se level of the IDAC system 
is monitored in each test by calculation of the equivalent temperature 
of an unused and shorted channel as shown in Figure 30* This level is 
always less than ■♦•O.OOSOC* 

3. Dynamic Operation 

In the transient thermal condition, control effectiveness must 
be defined by considering the influence of the container and sample 
mass thermal inertia. In the design of the present container, a 
trade-off was made which sacrificed the low thermal inertia of the 
container for the ability to determine the heats of phase change, the 
thermal conductivity and the specific heats of the sample mass* This 
was done by adding sensors in each wall of the container* The ability 
to obtain these properties makes the system unique and extremely 
attractive as a research tool for the development of protocols for 
freezing and thawing biological tissues* 

Typical real-time recorded plots of temperature and power versus 
time for distilled water are shown in Figures 33 and 34. These are 
the first-order calibrated data. The plots were made on an X-Y 
plotter with Y as temperature or power and X as time* The curve 
identified as (T 0 (i)) presents part one of the sample mass 
temperature-time history* At the end of this plot the recorder was 
reset so that the final time history could be plotted. This is 
labeled as (Tq( 2 )). The curve (T(,(]^)) presents the cooling 
history from +20^C to -40®C at ll^C/min and the curve (Tq( 2 )) 
presents the subsequent warming history, returning to the initial 
temperature at the same rate. In a similar manner the curves 
(Pa(i)) and (Pa( 2 )) present the initial and final histories of the 
power dissipation of the (A) face heater P^. The curves are 
qualitative and meant to present the overall view of the test. 

Curve (Tq(i)) shows an initial hold at a constant temperature 
of 20^C followed by a cooling ramp of 20^C/min. During the 
cooling ramp there is a sudden increase in temperature associated with 
a reduction in power. This is caused by the freezing of water* an 
exothermic process- Note from the position of the beginning and 
ending points of the ramp that the power curves (P^d)) ^rid 
(Pa( 2 )) shifted one unit to the right of the temperature 



curves « Consequently the deviations of temperature and power occur j 

simultaneously. This is an indication of the excellent response | 

between the heaters and the sensor within the sample mass even though i 

direct control of the heaters is by the sample mass surface j 

temperature as indicated on the curve. Upon resetting ! 

the time axis of the recorder aft®r the temperature of -40^0 is I 

atwiained and held for 120 seconds » the warming history of the sample i 

mass center is presented as curve (Tq( 2 )). Once again the change of 

phase is noted for both the temperature and power curves. This occurs | 

during the melting process which is endothermic, therefore requiring 
additional heat to maintain a steady temperature rise* Finally the 
temperature is held at the starting value of +20®C. 

Figure 34 presents the temperatures on either side of the sample 
mass for the same experiment presented in Figure 33, These are the 
temperatures under the direct cortrol of the heaters. The temperature 
and power curves are shifted one horizontal unit from each other on 
the graph. Note the excellent linearity of the temperature ramp 
during cooling, and the elimination of the phase change temperature 
deviation from this ramp during the warming process. In the freezing 
process a small control loss becomes apparent because of the 
suddenness of the phase change from the supercooled state. In some 
tests where supercooling is minimal the freezing perturbation as seen 
here is not evident from the curve. Some investigators have claimed 
that the temperature of a finite mass can be controlled during the 
phase change. As is shown in Figures 33 and 34, this may be apparent 
for a portion of the mass, e,g,, the surface in this case, but because 
of the finite thermal conductivity of the mass, it can not be true in 
general , 

Figure 35 presents a thermodynamic cycle of the same test given 
in Figures 33 and 34. The temperature at the center of the sample 
mass (Tq) is plotted versus the average power of the heaters. Point 
one is the starting point of the cycle; the temperature is held at 
*f20Oc. As power is decreased the temperature is reduced in a 
non-linear manner until the thermal inertia of the system and mass is 
overcome. When this has occurred, there is a linear reduction of 
temperature with power until freezing of the sample mass begins. The 
actual freezing point of distilled water is 0®C, however 
supercooling can delay freezing until about -10®C, At this point 
there is a rapid decrease in power with an associated increase in 
temperature. As the phase change reaches completion, there is an 
increase in power to drive the sample mass temperature back to the 
controlled ramp followed by a second nearly linear ramp. Point two is 
the second hold point in the test and is obtained by an increase in 
power with little change in temperature again associated with the 
system thermal inertia. After holding for two minutes at point 2, the 
warming cycle was started. There is a nearly linear increase in 
temperature with power until about -5®C where thawing becomes 
evident and an increased power is required to provide the additional 
heat for melting. Finally, the system returns to the starting point 
and the test is concluded. 
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The area within the primary envelope of the cycle presented in 
Figure 35 is the hysteresis loss due to thermal inertia of the sastple 
mass and its container. This area is a function of the starting and 
ending temperatures and, for a given sample mass, the rates of coding 
and warming. This loss is described by measuring the lag between the 
control system's reference temperature and the response of the point 
in the sample mass of least sensitivity to temperature control, i.e,, 
the sample mass center. To describe this dimension of the envelope 
the thermal inertia power (TIP) is defined as the average power 
measured at constant temperature for a given cooling or warming 
protocol. This measurement disregards deviations of the envelope due 
to phase changes. Figure 36 presents the relationship of TIP as a 
function of cooling or warming rate for distilled water. At present, 
the only significance for the TIP is that of a performance 
characteristic of the entire system. For detailed analysis, however, 
it may be a significant parameter in dealing with the sample mass 
performance during phase change. The relationship of the TIP to the 
temperature rate appears, within the accuracy of existing data, to be 
the same for all sample masses considered thus far. This includes 
distilled water, 0«9Z saline solutions, and packed red blood cells. 

Figure 37 shows a method for describing the performance of the 
control cycle for the same test discussed in Figures 33 to 36. It 
describes a "worst case" condition, since it presents the difference 
between the reference temperature of the controller and the 
temperature of the center of the sample mass. Deviations from the 
primary envelope are due to phenomena resulting in maximum loss of 
control, i.e., phase changes. The general shape of the primary 
envelope is an indicator of the thermal inertia of the system. 

Because the data in the figure are not corrected as shown in Figure 
31b, the value of the function does not remain at zero but varies to 
positive and negative values. This plot is interesting because it 
tends to amplify the phase change as an event in the cycle and may be 
useful for the determination of other events during the detailed 
analysis of red blood cells or other biological media. 


F. WATER TEST RESULTS 

1. Freezing Water 

The freezing of water is an exothermic (heat release) process 
that theoretically occurs slightly below the ice/water equilibrium 
temperature of 0®C. Pure bulk water, in the absence of ice crystals 
however has to be supercooled to -lO^C before freezing occurs. 
Ordinary water containing foreign bodies freezes from about -6^C to 
-25<^C. 


Once ice begins to form, the latent heat must be removed from 
the sample to insure continued phase '^hange. This is the cause of the 
temperature excursions seen in Figures 33, 34, 35, and 37. This heat 
removal may occur from heat trau‘^fer a) through the supercooled liquid 
phase under a temperature gradient directed toward the freezing 
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interface or b) through the solid phase under a temperature gradient 
in the ice. Any motion of the ice/water front is conditioned by 
thermal conduction of the system to a heat sink. Cooling rate and 
rate of ice formation is determined by the temperature gradient in Che 
water or the ice. Both cases mentioned above are seen to occur in the 
experimente presented in the following figures. 


2. Test Data For Distilled Water 

Freeaing and thawing tests were conducted with pure distilled 
water to arrive at a better understanding of Che behavior of the 
system during temperature control and the sample mass during phase 
change. Figures 38a to 38d show data recorded during an experiment in 
which the temperature of the water was reduced in steps from -»-20^C 
to -60°C as noted on Figures 38a and 38b, and subsequently increased 
back to -f20^C at a constant rate of 4'2°C/min. as shown in Figures 
38c and 38d. Each step was held for a period of 2 minutes and cooling 
occurred between steps at a rate of li^C/minute. Two curves are 
plotted on Figure 38a. Tq is th ' temperature at Che center of the 
sample mass (water) and Tq (1), "q (2) and Tq (3) present a 
sequential temperature history ' >r this point. is the power 
dissipated in Che (A) side heater and is also presented sequentially. 
Each point on Pa(D corresponds in time to each point on TgCl) 
witlt the Pa curves shifted one unit to Che right. Since these 
curves are for qualitative study, the vertical axis is not calibrated, 
although temperature can be inferred from the hold points. In all 
cases, each major unit on the horizontal axis is 50 seconds as 
indicated in Figure 38a. 

Freezing occurred during the cooling between the steps at 0°C 
and -5°C. This is shown in Figure 38a by noting the increased 
temperature in the Tgd) curve and the decreased power in the 
Pa( 1) curve. Because of the tight temperature control of the Tax 
and Tgx sensors, this phase change is not apparent in Figure 38b. 

This was not the case for the experiment described in Figure 34 where 
both Tax and Tgx indicated phase change at the same cooling rate 
for the curve of Figure 38a. This difference is attributed to the 
fact that the Figure 34 test allowed the water to supercool to -10°C 
before freezing. The associated change in the sample mass temperature 
was too sudden to allow control. This is seen by the sudden drop in 
heater power as shown in Figure 33. In contrast, the apparent 
temperature control of Tax ‘^Bl indicated in Figure 38b and 
the small deviation in the value of T^d) in Figure 38a resulted 
from less of a demand on the heater as indicated by the more gradual 
changes in the heater power PA(i) sho%m during the phase change 
in Figure 38a. Consequently the intuitively predictable result is 
that the lesser the degree of supercooling prior to phase change, the 
greater is the apparent temperature control of the sample mass. The 
term "apparent temperature control" is used because a temperature 
gradient must always exist in the sample mass in order to exchange the 
heat of phase change. Nevertheless, allowing this change of phase to 
occur over a longer period of time as was the case shown in Figure 
38a, appears to reduce the temperature gradients within the sample 
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mass and conaequently allows Che appearance of being under better 
theraal control. Because of published evidence that indicates reduced 
freezing damage to cells undergoing phase change under minimal 
supercooling, the control and measurement of the supercooling/phase 
change phenomena as just discussed is of considerable importance in 
the study of cellular freezing. 

Figure 38e presents a plot of Che center temperature, To, 
versus the average heater power. Points 2-11 indicate points where a 
constant temperature/ power setting was held for a period of 120 sec as 
indicated in Figure 38a. Point 1 is the start of the test and point 
12 is Che end. The hold points do not fall on the cooling curve due 
to Che thermal mass of the container and the sample (water). The 
freezing process occurs between points 3 and 4. Thawing occurs in 
nearly Che same temperature range on the warming curve. 

It has already been mentioned that the test of Figures 38a-38d 
showed freezing under a reduced supercooling condition compared to 
that of Figure 33. This is shown more clearly by considering the same 
two tests as indicated by Che comparison of Figure 35 with 38c and 
Figure 37 with 38f. It is readily apparent from Che temperature 
excursions during freezing indicated in Figures 38e and 38f, that this 
case experienced a reduced gradient within the sample mass. This 
condition is easier to note by comparison of the temperaCure/power 
cycles Chan by the temperature or power time histories, and presents a 
valid agrument for consideration of the cyclic process plots. 

The remaining water test experimental data are presented in 
Appendix D. There is no discussion of this data but they are included 
for Che purpose of completeness of this report. They add no new 
concepts to those already presented but they do support the 
conclusions obtained on the performance of the system. 


G. SALINE TEST RESULTS 

1. Freezing Saline 

VHien ice forms in a salt solution, Che crystals are pure water. 
As the temperature is lowered below this, the dissolved salts become 
concentrated until all the water that can be crystallized as ice 
freezes leaving only the solutes and their water of hydration. 

Further reduction in temperature '.-"suits on the solidification of the 
remaining solution. The lowest t"^v®*‘«fure at which the solution 
remains liquid is called the eutectic temperature. For NaCl 
solutions, this is -21.8°C. The concentration of solute necessary 
to achieve this minimum temperature is the eutectic concentration. 


2. Test Data For Saline 

Figures 39a through 39f present data typical of the freeze/thaw 
tests conducted with 0.9% NaCl solution (physiological saline). As in 
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th« dltcuaiioo of th« distilltd water axpariaMnCit tUa data ara 
praaantad in plots of actual teat data and theraodynamic cycle of 
taaparatura varsua power. 

Piguraa 39a to 39d present the tine hiatoriea of Tg, 
and Tgi for the cooling and heating phases of the test. The 
test is sinilar to that accomplished for distilled water with the 
exception that the -60^ hold point is not included in the plots. 

The hol^ points were for a duration of 120 sec. and the cooling rate 
between hold points was ll^/nin. As in the water experiments y the 
warming rate was at 2<>C/min. 

In Figure 39a the exothermic freesing process is seen to occur 
during the -10<^ hold point. The sample mass center temperature 
(Tg) rise during Tg(2)y occurs simultaneously with the sudden 
power decrease shotm during Pa( 2). After this occurrenesy there is 
no further indication of sudden phase change activity. 

Figure 39b shows the and Tg]^ temperatures during the 
cooling phase of the experiment. The phase change was not app- ■> - on 
this plot. It is obvious from Figure 39a however that the samp..: lass 
temperature was not held at a constant -lO^C as isay be erroneously 
concluded from Figure 39d but varied considerably as indicated b** the 
gradient between Tq of Figure 39a and Tgx or Tgx of Figure 39b. 

Figure 39c presents the Tg and time histories for the 

warning phase of the experiment. The endothermic phase change occurs 
slowly over a period of about 200 sec. During this time the 
temperature is seen to drop below the desired control ramp indicated 
by the dashed liney and the power is seen to slightly increase above 
its nominal ramp for the given warming rate. The temperature decrease 
below the desired ramp is a laeasure of the gradient across the sample 
mass since no apparent temperature deviation appears at the control 
sensors T^i end Tgx ** seen in Figure 39d. 

Figure 39e presents the sample center temperature Tg as a 
function of the average heater power in a fashion similar to Figure 
38c. The freexe phase change occurs at the -10°C hold point and the 
thaw appears <is a subtle bend in the warming curve near O^C. 

Figure 39f i^ analogous to Figure 38 f, representing the 

temperature differential bet%reen the sample mass center temperature 
(Tg) and the reference temperature which the system is attempting to 
reach. The large difference in this value during the non phase change 
process is because the data calibration is approximate and the 
corrective equation (Figure 31b) was not applied. The important point 
of the plot is the large excursions which occur during the freeze and 
thaw phase changes, once again indicating that this correlation may be 
an excellent method to analyze the phase changes of a complex 
biological media. 
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The remainder of the saline test data are presented in Appendix E 


H* RED BLOOD CELL TEST RESULTS 

1« Freezing Red Blood Cells 

The quality of freeze-preserved blood products has improved 
significantly over the past 10 years* The increasing use of these 
products in medicine has been encouraging but still falls short of 
their full potential. The high cost of processing red blood cells 
through the freeze preservation protocol and with the use of 
cryoprotectants is most likely the reason for the lack of use of this 
method* In an attempt to reduce this cost^ a system has been built 
for the conduct of research on red blood cells without the use of 
cryoprotectives * 

The literature on the freezing of red blood cells is quite 
complex and, in some cases, contradictory* Consequently, no attempt 
is made in this study to describe the biochemical and mechanical 
changes associated with the cells but to show that an effective device 
has been produced to accurately control the temperature-time history 
of the cells and to yield certain thermodynamic properties of the 
cells as a function of temperature* 

2. Test Data For Red Blood Cells 

Fresh blood was drawn in 10 cc vaccutainers containing 1 cc of 
acid citrate dextrose (ACD)* The blood was centrifuged in the 
vaccutainers at 3000 rpm for 3 minutes to separate the red cells. 

Cells to be used in each test were aspirated from the center 5 cc of 
each vaccutainer. The cells were transferred to the test apparatus by 
syringe. The container and transfer apparatus were flushed with 0.9% 
saline solution and finally with ACD prior to cell transfer* 

Potential cell damage due to the transfer of packed cells with 
this method was studied by aspirating the cells from the 
freezer/container after allowing the cells to rest in the container 
for 30 minutes at room temperature. After aspiration* the cells were 
placed in a vaccutainer rinsed with ACD and containing 5 cc of 0*9% 
saline* After mild agitation, the sample was centrifuged at 3000 rpm 
for 3 minutes* No noticeable hemolysis was indicated, thereby 
assuring minimal cell meclianical damage due to the transfer protocol* 

Figure 40a presents the center temperature (Tq) and the A 
surface power (P^) histories of a typical test with packed red blood 
cells. The initial temperature hold point was set at 15®C. A 
cooling ramp at lO^C/min. was conducted until a temperature of -5® 
was net. The temperature of -5^C was held for nearly 350 sec* 
during which the phase change occurred* A second cooling ramp at 
lO^C/min. was then conducted to a temperature of -40^C* This 
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temperature was held for 120 sec. after which a warmirg ramp df 
lO^C/min. was conducted to a temperature of -l^C. Prior to 
reaching this point the ice began to melt. This is seen by the slow 
change in the Tq (3) curve as it approached -1®C and by the slow 
return of the power curve (3) to a steady state condition. 

During the 400 sec. hold at -l^C, the phase change continued. 

Finally, a second warming ramp at 10^/min. was established. This 
stimulated the final melting as indicated by the nonlinear temperature 
curve and the increased heater power. The test concluded at the 
temperature of ■•■15^0. 

The response of the sample volume center temperature to control 
and the obvious occurrences of phase change during freete and thaw are 
excellent examples of the application of the control system to the 
thermodynamic study of freeze preservation of biologicat media. 

Figure 40b presents the sample volume center temperature versus 
the heater average power in a thermodynamic cycle for the entire test 
presented in Figure 40a. Mote that the phase change occurs at the 
hold point (2). Cooling then continued to the second hold point at 
(3). The thawing of ice begins to appear during the warming phase 
prior to reaching the hold point (4) as sho%m by the nonlinear 
behavior of the curve. After leaving hold point (4), the non-linear 
behavior continues to show a phase change. The test ends at the 
starting point, (1). 

Figure 40c is a thermodynamic cycle for the same test but 
presenting the difference between the sample center temperature (Tq) 
and the controller reference temperature (Tggp) as a function of the 
average heater power. This differential temperature is an indication 
of the dynamic tempei'ature control of the system and is exaggerated 
during any phase change. Again note the excursions which occur during 
the hold at (2), and immediately before and after the hold at (4). In 
this case the excursion prior to the hold at (4) is even more obvious 
than the preceding graphs. This cycle analysis approach as applied to 
the identification and location of phase changes, may prove to be 
quite useful for future research in freeze preservation techniques. 

Similar plots for the remainder of the red blood cell tests are 
presented in Appendix F. 

Approximate tests for the presence of hemolysis were conducted 
on some of the experiments. These were accomplished by mixing the 
aspirate from the test apparatus with an equal volume of 0.9Z saline, 
centrifugation at 3000 rpm for 3 min., and visually checking the 
supernatant for discoloration. The results were categorized into 3 
qualitative descriptions of severe, moderate and slight as indicated 
in Figure 41 . The intent of these tests was not to prove or disprove 
the efficacy of a particular protocol but to show differentiation for 
the tests where freezing damage was obviously present. Note from 
Figure 41, cases 7, 13, 14 and 15, that the mere occurrance of a phase 
change was not sufficient to result in severe hemolysis. In cases 
where no phase change occurred, however, as in 1-6, 8, and 11, slight 
or no hemolysis was always indicated. 
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SECTION III 




CONCLUSIONS 


Considerable effort has been made to understand the performance 
of the system* This has been done in order to produce a precise 
method for the study of freeze/ thaw protocols of biological tissues* 
These tissues exhibit complex and in many cases ^ unknown changes with 
temperature. The ability to make accurate and precise measurements of 
the heats of phase changes as well as of the thermal conductivities 
and specific heats will relate directly to the understanding of these 
changes. The system that presently exists is not just a blood 
freezer. It is an instrument to study the cryogenic behavior of 
complex substances of organic or biological nature. 
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Figure 1. Basic Design Concept 
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Figure 2. System Block Diagram 
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PORT 

b = BOND 

h = HYDROPHOBIC GREASE 
= SPRING FORCE 

V = VESPELSP21 

R, _5 = ALUMINUM PLATES, HARD ANODIZED 

S, _3 = NICKEL TEMPERATURE SENSORS WITH KAPTON BACKING 
H = NICHROME HEATER WITH KAPTON BACKING 

X = LOCATION OF THERMOCOUPLE AT CENTER OF MASS TO BE FROZEN 
Sj = CONTROL SENSpRS 
S2_3 = HEAT FLUX SENSORS 


Figure 6. Model III Configuration 
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Figure 9. Progrannnable Freeze/Thaw Container for Biological Cells: 
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Figure 11. Programmable Frecze/Thav Container for Biological Cells 
View of Exterior Wall 
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Figure 14. Data System Temperature Control Block Diag 
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Figure 16. Temperature Controller, View 
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Figure 22. Card #4 D/A Converter and Sensor Excitation Supply 
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IDAC Functional Block Diagram 
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Figure 33. and Versus Time for Distilled Water (Typical) 
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Figure 36. Tlierraal Interla Power Versus Temperature Kate of Change 
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Fii^urc* 38c. T and P, Versus Time for Distilled Water Dtjring Warming ^Typical) 
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Figure 38f. to Refort^nce Temperature Differential Versus Average 

Power for Distilled Water (Tvplcal) 
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Figure 39b. T and T . Versus Time for Saline During Cooling (Typical) 
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Figure AOa. and Versus Time for Red Blood Cell Cooling and 

Warming (Typical) 
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Figure 40c. to Reference "emperature Differential Versus Average 

Power for Red Blood Cells (Typical) 



0 Phase Chianj?c Hecolysit. 



rrjnary of Pod Blood Coll Test Pesult; 



APPENDIX A 
PATENT DESCRIPTION 


A- 1 


1”^ ' 



1 

1 


H . , . ■■ Cl 


! 

\ 

\ 

i 

1 


U. Govw.ni'ierft; 



United States Patent ini 


(HI 3,733,463 


Low et al. 


|45| May 15, 1973 



134) TEMPFRATtRE CONTROL SYSTEM 
WITH A PULSE WIDTH MODULATED 
BRIDGE 

(76) Inventors George M. Low, Acting Adininistri- 
tor of the National Aeronautics and 
Space Administration with respect 
to an invention of. Richard C. 
Heyser, Tujunga. Calif. 

[221 Filed Dec. 24, 1970 

121) AppI No.: 101,214 

(52) U.S. Cl 219/490. 219/50 

[511 Inc. Cl H05b L02 

(58) Field of Search 219/404. 497. 499. 

219/501. 317/42. 323/75 H. 75 K 

[56) References Cited 

UNITFD STATES P.ATENTS 
2ol>.372 1/1958 Boc'kcr et al 219/499 


3.526.272 9/1970 Walts el al 2 ! 9/497 \ 

2.681.431 6/1954 Wannamakcr. Jr K 

3,159.796 12/1964 Van Sand>v>k 323/75 K 

3.546.434 12/1970 Apcl 219/499 UX 


Primary f^uminer— Bernard A. Gilheany 
AsMStant tiuminer— F F. Bell 

.^/fornev- Monte F. Mott, Paul F. McC'aul, John R. 
Manning 

[571 ABSTRACT 

A temperature control system is disclosed which in< 
eludes a modified wheatstonc bridge with a resistive- 
capacitive (RC) circuit in one leg of the bridge The 
RC circuit includes a resisior which provides an effer- 
tivc resistance as a function of its absolute resistance 
and the on-iime to off time ratio cf pulses supplied to 
a switch connected thcreacross A sawtooth voltage is 
priHluced across the RC circuit, the voltage being 
compared with the voltage across a tcmperatuie sen- 
sor. with heal being applied during each pulse period 
portion when the sawtooth voltage exceeds the voltage 
across the lemperaiuie sensor. 

8 CUims, 2 Drawing Figures 
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TEMPERATURE CONTROL SYSTEM WITH A 
PULSE WIDTH MODULATED BRIDGE 

ORIGIN OF THE INVEHnON 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public L^iw 
85-568 (72 Stat. 435. 42 USC 2457). 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention generally relates to a tempera- 
ture regulator and more particularly, to a s>‘stem for ac- 
curately controlling the temperature of an environment 
within an oven or the like. 

2. Description of the Pnor .Art 

A ty'pical temperature control svstem such as that de- 
scribed in U S. Pat. No. 3,136,877 consists of a conven- 
tional Wheatstone biidgc circuit. The bndge circuit in- 70 
eludes a temperature seruitive resistor located in the 
environment, whose temperature is to be controlled. 
When the temperature is too low, a first bndge unbal- 
ance condition anses which is detected by a sensor 
which closes a switch, thereby enabling power to be 75 
supplied to a heater which heats up the environment. 

The environment is heated up until the temperature ex- 
ceeds a desired nominal temperature. .As a result, a sec- 
ond bndge unbalance condition anses. causing the sen- 
sor to open the switch and thereby terminate the supply 

of power to the heater. 

With such a prior art svstem. the environment tem- 
perature can only be controlled within a finite sigmfi- 
canl range, since power supply to the heater occurs 
whenever the temperature drops below the nominal 
temperature and is terminated whenever the tempera- 
ture IS above the nominal temperature. Even wuh sensi- 
tive instruments, tnc temperature range is quite signifi- 
cant, and is often greater iluin that required for very 
precise studies or applications, such as in space explo- 
ration siudirs. Thus a need exists for an improved sy%- 
tern for accurate temperature control. 

OBJECTS AND SUMMARY OF INE INVENTION 

It is a primary object of the present invention to pro- 
vide a new improved temperature control s>*stcm. 

Another object of the present invention is to provide 
an improved temperature control system which incor- 
porates a novel wheatstone bndge circuit 

A further object of the present invention is lo provide " 

1 new system which is programmable to control the nse 
in and the maintenance of an environment tempera- 
ture. 

These and other objects of the present invenuon are 
achieved in one embodiment by incorporating in one 
leg of a wheatstone i^dge circuit fust and second sen- 
ally connected resistors, with the second resistor being 
connected across a switchuig transistor. The traniistor 
is successively switched between on and off states, dur- ^ 
ing which the second resistor is either shorted out or is 
in senes with the first resistor. During each switching 
cycle the effective resistance of the second resistor 
which is in scries with the first resistor depends on the 
ratio of the transistors on-time to off-time. By varying ^ 
this ratio, namely by pulse width modulating the 
switching of the uansistor, the effective resistance can 
be changed. 
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In acoordance with the teachings of the present in- 
venuon. a capacitor is connected across the two resis- 
tors. During each switching cycle the voltage across the 
capacitor which is in the form of a sawtooth is com- 
3 pared with the voltage acrxxs the temperature sensor in 
the bndge. During each cycle, power is applied to a 
heater only if the capacitor voltage exceeds the temper- 
ature sensor in the bndge. 

The novel features of the invention are set forth with 
10 particularity in the appended claims. The invenuon will 
be best understood from ihc following descnption 
when read in conjuncuon wi Ji the accompanying draw- 
ings. 

j 5 BRIEF DESCRIPnON OF THE DRAWINGS 

FIG. 1 is a viiagram of a basic embodiment of the 
present invention, and 

FIG. 2 IS a multilme diagram useful in explaining the 
embodiment shown in FIG. 1. 

DESCRlPnON OF THE PREFERRED 
EMBODIMENTTS 

Attention is now called to FIG. 1 which is a combina- 
tion block and schematic diagram of the present inven- 
tion. Therein, numeral 10 represents a wheatstone 
bndge arcuit, connected across a source of potential, 
such as a DC battery 12. The bndge 10 irKludes a resis- 
tor R1 which is connected in series with a temperature 
sensitive resistor or simply a temperature sensor Pt 
across the battery 12. It is sensor Pt which is located iri 
the environment or oven whose temperature is to be 
controlled, and whose resistance varies as a function of 
temperature. The junction point X between Rl and Pt 
is coruiected to one input of an operational amplifier 
IS. The bridge also includes another resistor R2 which 
is serially connected with variable resistors PI and P2 
across battery 12. A capacitor C is connected in paral- 
lel across PI and P2 and the junction point Y between 
R2 and PI and C is connected to tire other input of am- 
plifier 15. 

The latter is operated as a comparator to provide ^r. 
output at point Z of a preselected level which closes a 
switch 17 only when the voltage at point Y is greater 
than that at point X. Wlkcr. the switch 17, which is 
represented in FIG. 1 by a transistor Ql, is closed it 
connects an oven heater clement 20 across a source of 
potential, shown in RG. 1 at 4-28V. Thus, only when 
switch 17 is closed is heat supplied to the oven by ele- 
ment 20. 

As shown in F!G 1, connected across the resistor P2 
is a switch 22. in the form of a transistor Q2, which is 
activated by pulses such as pubes 26 and 27, from a 
source 25, which arc supplied at a selected rale. Each 
pulse, such as pulse 26, has an on- time 26u during 
which the switch is closed so that P2 is shorted out. and 
an off-time 266, during which the switch is open so that 
P2 is not affected. ITie ratio of on-time lo off-timc is as- 
sumed to he controlled by source 25, which may in turn 
he controlled by a programmed computer (not shown). 
The ccntr\>l of tl.e cn-time to off-time ratio may be 
thought of as controlling or imxJulanng the pulse width 
of each pulse. Therefore, source 23 can be thought of 
as a source of pulse-widih-modulatcd (PWM) pulsca. 

In effect, during the on-time of each pulse P2 b 
shotted out, thus its resistance is rtro, neglecting the 
collector to emitter resistance of 02. On the other 
hand, duniig pulse off-time the full resbiance of P2 is 
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in scries with PI. Thus, by controlling the on-time to 
ofT-time ratio the efTecuve resistsjice of P2 during the 
full pulse duration is controlled. For example, if the on- 
time is only 25 percent of the pulse period, the cfTective 
resistance is one-fourth the resistance of P2. On the 
other hand, if the on-ome is half the pulse period, the 
effective resistance is half the P2 resistance. In practice 
PI is selected to equal the resistance of Pt at an ambi- 
ent temperature, while the effective resistance of P2 is 
chosen to equal the change (increase) in the resistance 10 
of Pt due to a change from the ambient temperatui€ to 
a chosen temperature to be controlled. 

The voluge at point Y changes due to the change in 
the resistance which P2 provides in series with PI. In 
the absence of C the voltage at Y would have a square 15 
wave shape as shown in line a of RG. 2. The upper volt- 
age level 31 would be present when Q2 is off and the 
lower voltage level 32 would be present when Q2 is on 
and P2 is shorted out. Therein, it is assumed that the 
off-time to on-time ratio is about 1;3. That is, during 20 
each pulse penod the transistor Q2 is on for about 75 
percent of the period. However, with C in the circuit 
due to the time constant of the circuit formed by C, PI 
and P2, the voltage at point Y during the period of each 
pulse from source 25 does not have the square wave 25 
shape shown in line a but rather the sawtooth shape as 
shown in RG. 2, line b and designated Vj. 

From the foregoing description it should thus be ap- 
preciated that in the present invention, during each 
pulse period, heat is applied to dement 20 only during ^ 
that portion of the pulse period when the voltage at 
point X designated in RG. 2 by the dashed line is 
less than Vy. Assuming that is less than \y during 
an entire pulse period, power is supplied to heater ele- 
ment 20 during the entire period, as represented by nu- 
meral 35 in line c which is used to diagram the power 
profile supplied to heater element 20, and which corre- 
sponds to the output of amplifier 15. On the other 
hand, if is greater than the maximum value of Vy 
during an enure pulse period, the heater element is ^ 
completely off during the entire period. Thus, the 
amount or quantity of heat which is supplied by heater 
element 20 to the oven during each pulse period is a 
function of the actual oven temperature as represented 
by Vy and the shape of Vy which is controllable by the ^ 
time constant and the pulse width modulation. 

From uic forcgcir.g it «hnuld be appreciated that, as- 
suming a constant pulse frequency, by conuuML'.g 
pulse width modulation, the rate at which the oven 
temperature is raised to a desired temperature may be 
controlled. Furthermore, once this ten ~>erature is 
reached it may be accurately maintained by controlling 
the pulse width modulation so that during each pulse 
period only that amount of heat which is necessary to 
maintain the oven at the desired temperature is pro- 
vided by heater rlcment 20. Maximum *ate of tempera- 
ture rise may be achieved by making the on-time per 
pulse penod equal to zero thereby raising the level of 
Vy to a maximum so that power is supplied to the ^ 
heater element dunng the full period of each pulse. The 
rate at which temperature is increased from a first tem- 
perature such as rl to r2 may be controlled by the pulse 
width modulation of the pulses to produce a desired 
rate at which power is supplied to the heater element ^ 
20. Once the desired temperature is reached, the nul«e 
wtdUi modulation is selected to insure the supply of 
power per pulse period which is required to maintain 
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the oven at the desired temperature. Experience with 
the novel invention has demonstrated the ability to 
control an oven temperature to within less than 0.01^, 

e g., o.ooix:. 

Referring again to RG. 2, let it be aanunod that the 
final desired temperature is that represented by the line 
36 for VjT and that the power or h^ per pulse needed 
to maintain that temperature is represented by numeral 

37. Let it fiinhcr be assumed that the oven is below the 
desired temperature. Thus, more power than is re- 
quired to maintain the desired temperature is supniied 
to the heater per pulse, as represented by lines 35 and 

38. In practice, the oven temperature may briefly oscil- 
late above and below the desired temperature until the 
temperature is stabilized. In RG. 2, lines 39 and 40 rep- 
resent Vj when the temperature is above the desired 
temperature. In response to the level represented by 
line 39, less than the required heat is supplied as repre- 
sented by line 41. Since level 40 is above the maximum 
value of Vy no power is applied to heater element 20 
as represented by line 42. Line 44 represents an oven 
temperature below the desired temperature. Conse- 
quently, more power than necessary to maintain the de 
sired temperature is applied as represented by line 45. 
L ines 48 and 49 respectively, represent at the de- 
sired temperature and the power per pulse period 
needed to maintain the desired temperature. 

Although particular embodiments of the invention 
have been desenbed and illustrated herein, it is recog- 
nized that modifications and variations may readily 
v'^cur to those skilled in the an and. consequently, it is 
imonded that the claims be interpreted to cover such 
mocLflcations and equivalents. 

Wh/t is claimed «s: 

1. A temperature control s>*stem comprising: 

a wheatstone-like bridge having fl.st and second par- 
allel legs and including a temperature sensitive ele- 
ment, whose resistance varies as a function of tem- 
perature, connected in series with a first bridge re- 
sistor in said first leg, and a resistor-capacitor cir- 
cuit connected in series with a second bridge resis- 
tor in said second leg; 

a source of power, 

heating means for providing heat when connecte i to 
said source; 

means for applying a voltage across said bridge, 
whereby the voltages across said element circuit 
are functions of the resistances therear^ross, respec- 
tively, 

means coupled to said elemcru ar.d to circuit for 
connecting said heating means to said source when 
the voltage across said circuit is greater than the 
voltage across said element; and 

control means coupled to said circuit for controlling 
the resistance across said circuit to vary in a prese- 
lected m«anner during each of a succession of pulse 
penods with the capacitor of said circuit control- 
ling the voltage of said circuit to vary in a prese- 
lected manner from a first level to a second level 
during a first portion of each pulse period and from 
said second level to said first level during the rest 
of the pulse period following said fust portion. 

2. The arrangement as recited in claim 1 wherein said 
circuit consists of first and second resistors connected 
in series and said capacitor is connected in parallel 
acron at least one of said resistors and said control 
means is coupled to said secona resistor and includes 
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first means for shorting out said second resistor during arrangefnent compti i ing: 


Mud first portion of each pulse period, with said first 
portion being selectively variable. 

3. The arrangement as recited in claim 2 wherein said 
control means include second means for providing a s 
succession of equal period pulses, each pulse being of 

a first level during said first portion of each pulse pe- 
riod and of a second level during the rest of the pulse 
period, and means for applying said pulses to said first 
means to short out said second resistor during the first io 
portion of each pulse period. 

4. The arrangement as recited in claim 2 wherein said 
capacitor is connected in parallel across both said first 
and second resistors, whereby when said second resis- 
tor is shoned out the voltage across said capacitor, 15 
representing the voltage across said circuit, rises gradu- 
ally from said first level to said second level and falls 
gradually from said second level to said first level dur- 
ing the rest of each pulse period when said second resis- 
tor is net shorted out. with the voltage across the aicuit 20 
during each pulse period haiing the shape of a saw- 
tooth. 

5. The arrangement as recited in claim 4 wherein said 

control means include second means for providing a 
succession of equal penod pulses, each pulse being of 15 
a first level during the first portion of each pulse period 
and of a second level during the rest of the pulse pe- 
nod, and means for applying said pulses to said first 
means to short out said second resistor during the first 
portion of each pulse period. 30 

6 . In a temperature control system of the rype com- 
prising a bridge with a temperature sensitive elcmcri 
whose resistance varies as a function of temperature in 
one leg of sauJ bridge, said s)‘sieni turther fociud«ng 
means for comparing the voltage acr.na said element 35 
with the voltage across another leg ef said brvige, an 


a remtor having flm and second terminals in said 
other leg; 

control means for oontrolling the resisunce acrou 
said terminals to vary from the full resistance of 
said reaistor to substantially zero during a selected 
portion of each of a accession oC periods, whereby 
the effective resstance across said terminals during 
each period is less than the resistance of said resis- 
tor, and 

a capacitor in said other leg connected in parallel 
across at least tiid resistor whereby the volLige 
across said other leg iricreases toward a peak value 
when the resistance across said terminals it the full 
resistance of said resistor an*i decreases from said 
peak value when the resistance across said termi- 
nals ui substantially zero. 

7. The arrangement as recited in claim 6 wherein said 
control means includes a switch which U connected 
across said first and second terminals, said switch being 
twitchable to an on state in which said terminals arc ef- 
fectively shoned out and an otT state during vihich the 
resistance across said tcmaiak is the resistance of said 
resistor, and source means coupled to eaid switch 
means for camtxolliiig said switch *0 be in the on state 
duiing a selected portion of each of said periods and m 
the off state during the rest of each of said periods. 

8 . The arrange n-.ent os recited in claim 7 wherein said 
source means is a source of a succession of equal pe- 
nod pulses, each pulse being of a first level which 
.twitches said switch U> iu on sute during a vanably se- 
lecied penod portion and of a second level which 
switches said switch to off state during the rest of the 
pulse txnod. 
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APPENDIX B 

CENTRIFUGAL RED BLOOD CELL FREEZE/THAW CONCEPT 
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In the initial conceptual design phase of this work a parallel 
design effort was established which considered both the flat plate 
container as discussed in the bulk of this text and the centrifugal 
device which is described in this appendix. 

The centrifugal system can be divided into three groups of 
components as shown in Figures B-1, B-2 and B-3. Figure B-1 depicts 
the gaseous nitrogen (GN 2 ) temperature control system. Liquid 
nitrogen (LN 2 ) gaseous nitrogen (GN 2 ) is passed through an 
electrical heater which establishes an environment temperature as 
measured by the thermocouple rakes to within about l^C of the 
desired blood mass temperature. The mixers are a set of baffles which 
produce an isothermal flow cross section. The flow then passes 
through a flow laminator to provide uniform convective heating or 
cooling of the cylindrical rotating tube which contains the red blood 
cells. The red blood cell mass is divided into a test volume, 
accompanied by adjacent guard volumes, so that maximum control and 
information of the test volume dynamics can be obtained. The tube is 
rotated by a stepper motor to allow uniform circumferential heat 
exchange with the flow. Electrical connections between the electronic 
control and monitoring system, and the temperature sensors and heaters 
on the tube itself are made through the mercury (Hg) slip rings at 
either edge of the tube. 

The tube assembly which contains the red blood cells is depicted 
in Figure B-2. The test volume of cells is separated from the guard 
volumes by spacers which can move in the longitudinal direction to 
accommodate freeze expansion. The thickness of the blood annulus 
within the tube is less than the spacers; consequently, radial freeze 
expansion is also possible. The fine temperature control of the cell 
mass is accomplished by a pair of bifilar heater/sensor windings, 
separated by a thin Kapton film which acts as an electrical 
insulator. These are shown in detail in Figure B-3. 

Figure B-4 depicts the block diagram for the system. A 
temperature profile memory is fed to the sequencer comparator which in 
turn senses the temperature of the cell mass and controls the heater 
power so that the input temperature profile is followed by the cell 
mass. The sequencer comparator also sets the heat sink (or flow) 
temperature by adjusting the GN 2 /LN 2 niixture ratio and by setting 
the GN 2 heater power. 
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Figure B-1 . CN Temperature Control System 
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APPENDIX C 


THKR>L\L ANALYSIS OF RED BLOOD CELL FREEZING SYSTEM 

In order to provide an approximate thermal analysis of the freeze/ 
tlinw system, the component structural elements of the container were 
reorganized into the heat tr.ansfer analysis model shown in Figure C-1. 
Figure C-1 relates each element of the model with that of the actual 
container as described in Figure 7. For the purposes of this analysis 
it is not as important to identify the real thermal properties of the 
device as it is to demonstrate the procedures by which the thermodynamic 
properties and characteristics of the sample mass can be evaluated. 

This is accomplished assuming one dimensional heat flow tlirough tlie 
container walls and neglecting the effects of radiation and convection. 

An energy balance is derived for each element of the analytical 
model except for the exterior elements. Tlie temperatures of these 
elements are used as boundary conditions for the analvsis and the ther- 
mal coupling ot these eli‘ments with the cold sink is uninformative for 
the purposes ol this work. 

When the energy balances are described for the general transient 
condition, tin* relationships for the thermal resistances between each 
of the model elements an* derived in terms of temperatures and powers 
monitored bv tlie existing system. Having these plivsical characteristics 
ot tile container, i cjuat ions a»'e developed for tlu* temperature dependent 
tliermal conductivity and specific *ieat of the sample mass, and for the 
sample mass hc-at of phase cliange 
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Figure Cl. heat Transfer Analysis Mod«*l 
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LIST OF SYMBOLS 



M 

X 

P 

X 

*^(X-Y) 

t 

T 

X 

T 

MN 

0 

(I 

0 


* 


area of element normal to heat transfer 
specific heat of element x 
correction factor (equation 41) 
correction factor (equation 45) 

instantaneous rate of heat production or absorption 
due to phase change 

thermal conductivity of element x 

mass of element x 

power generated in element x 

thermal resistance between elements X and Y 

thickness of element in direction of heat transfer 

temperature of element x 

average temperature of elements M and N 

t ime 

experiment function (equation 40) 
experiment function (equation 44) 

at a steady state condition with uniform temperatures 
between heaters 

at a steady state condition with an imposed temperature 
gradient between heaters 


1. Basic Knergy Balance 


The energy balance for a model element is formed by equating the 
rate of change of the internal energy of that element to the power 
dissipated within the element less the lieat conducted from the clement. 
In general tliis is written as; 
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g. Element (B3) 


Boundary Condition: Known as lunction of time 


2. SOLUTION FOR THERMAL RESISTANCE VALUES 


For a given steady state condition, the derivative of temperature 
with respect to time is zero, and the electrical power and temperature 
dil ferent ials are constant 


V = constant = 


AT = constant * AT* 


Using these values in Equation (1), there results the relation for 
the tliermal resistance as a function of temperature. 


\X-Y) p* 


if the resistance is not a strong function of t emi'erature , it can 
be approximated by noting that the temperature diiterential is the 
difterence of two known temperatures. 


AT = T -T. 
I J 


Tlierelore the approximate temperature level ol the element is 
calculated bv 


T = O.S(T. + I\> 
I 


The resistance is then expressed as a function of the approximate 
temperature level 







( 13 ) 


This basic approach, discussed above, is used to obtain the thermal 
resistances of the analysis model 


a. Thermal Resistance; R..,. 

iAZ-Ai ; 


Under steady state conditions, tlie temperature gradient between 

all elements internal to the heaters, i.e., T , T^, and T„, vanishes 

o A 1 lU , 

so that 
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o 




(U) 


Applying Equation (14) to the steady state form of Equation (3), 
y ields 


(T 
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where 
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(ib) 


Equation (13) states that the Mieimal resistance between elements (Al) 
and (A2) is an approxim.ite t unction of the average temperature ol those 
elements, l\ . . 

i\A 
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b. TlK'rmal Rosistanco: R 

(BJ-Bl) 

Kqualion (14) is applied to the steady state form of Kquation (4) 
to yield 


k (T *) = 

U52-B1)^ BA r * 


(17^ 


where 




(18) 


The elements of equations (17) and (18) are defined in a way similar 
to those of eijuations (15) and (In)* With the application in this case 
beiny; to the "B*' wall. 


c. Thermal Kesistance: K 


(A2-A3) 


The steady state form of luiuation (5) yields 


R /r * T *\ = R (T *) 

(A2-A3) \ AB ’ M J (A2-A1) ^ \A ’ 0\ -1\.J* 
' ' Ai AZ 


(19) 


where 
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AB 
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liquation (19) ^ives the resistance between elements (A2) and (A3) in 

terms ol tliat between (Al) and (A2) . Conseiiuent ly the dependent 

variable is a function of both the averay;e temperature between (A2) and 

(A3), given as T , iind that for the (Al) and (A2) elements, given as 
AB 

V . . . The presence of this unusual dependence is not a real physical 

i\A 
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dependence but simply results from the analysis approach. Inserting 
Equation (15) into Equation (UM there results. 
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Thermal Resistance: 


*^VB2-B3^ 


Using a method analogous to the development oi Equation (21), 
there results 


*^(B2-B3) 




^^B3 ^B2^ 

•'b* 


where 


BB 




( 22 ) 


(23) 


e. Thermal Resistance: R (tliermal conductivitv ol sample) 

(o-Al) 

The thennal resistance between the sample volume (o) and the 
interior wall (Al) is composed oi two elements. One element corresponds 
to the sample volume, the other to the interior wall 


(o-Al) 
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Al 


(24) 


In general an arbitrary therm il resistance is inverst'ly pi tq>or t i ona 1 ti' 
the therimil conductivity and the cross sect ional area oi the material 
ior whicii it is being applied and directly |M*oport ional to the 
thickness ol the material in the direction oi the heat Mow. 
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Then Equation (24) becomes 
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(o-Al) “ R A 
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For the device in question 



and the resistance becomes 
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Consi Jor ing tliat a Lvpical sample mass will have the conduct ivitv (K ) 

o 

oi watei" and Llie inner wall is aluminum (K ) , tl>en it is true tliat 
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SO that Kquation (,28) can be rewritten totally in terms of the 
sample mass 


K 


t 

o 


f o—A 1 ) K A 
o o 


K 


(30) 


Similarly lor the UO side 
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During periods when there is no phase change 


h - 0 
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( 32 ) 


Now a unique capability of the system will be involved. A steady 
state gradient is established such that 
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and 
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If the temperature difference across the sample mass is of the 

t 

order of a few degrees centigrade, then is a good representation 
of the sample mass temperature* Applying Equations (30), (31), and 
(32) to Equation (3) under the conditions of the steady state gradient* 
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Tlie thermal resistance ^(^^2-Al)* obtained from Equation (15) 
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Solving for the thermal conductivity of the sample mass from 
Equation (30) 
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Incorporating Equation (36) into Equation (37) 
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A similar equation can be derived tor the (B) wall. If this is done, 

k is determined as the average of the two values* 
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( 39 ) 


Equation (39) presents a valid theoretical prediction for the thermal 
conductivity of the sample mass. However, it often is advisable to 
reference the calculation to a mass of known conductivity such as 
distilled water. Given the same apparatus and the same value for 
(T^*), Equation (39) can be used to calculate the conductivity for 
water. This will involve the relationship given in Equation (39) plus 
a correction factor assumed to be a function of (T^^*). If Equation 
(39) is written as 


k (T*) = 
o o 



(40) 


where 4^ is the experiment function in the brackets of equation (39), 
o 

then the corrected equation for water is 


k (T*) 
w w 



(41) 
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where f(T*) is the correction correlation which also accounts for the 
w 

presence of Vespel in the sample volume. Applying the correction 
factor to the calculation of Equation (40) yields 


k (T') 

o'* o^ CORRECTED 




where (T*> 

w 


k (T') 

o' o*^ CORRECTED 




Equation (42) presepts the proposed algorithm to be used for the 
calculation of the temperature dependent thermal conductivity of the 
sample mass. 

3. SPECIFIC HEAT 

Having calculated the thermal conductivity of the sample mass, the 

value of R (equal to R, and R, _,.) can be obtained from Equa- 

o (o-Al) (o-Bl) 

cions (30) and (31). For experimental regions without phase change. 
Equation (1) can be solved for the thermal mass 


M C (T ) 
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The application of Equation (43) is theoretically correct as stated. 
However the normalization of the relation to water allows the applica- 
tion of a correction factor which is assumed a function of (T ). 

o 

Rewriting Equation (43) as 


M C (T ) - 4i (T ) 
o o o ^o o 
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Applying Equation (44) to measure tlie known specific heat of water, the 
correction factor B('l'y) is applied as with the calculation of the thermal 
conductivity. 
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(45) 


Solving for the correction factor and applying this to Equation (44) there 

results (letting T « T ) 
o w 


v\(T ) 

M C (T ) = M C ( T ) t-Tp'T 
OO O WW O (f(I) 

w o 


(46) 


Equation (4o) presents the recommended formula for the calculation of 
the specific heat of the sample mass, 

HEAT OE PIUSE ClL-VNtiE 


Having calculated the thermal conductivity and specific heat of 
the sample mass outside the regions of phase change, it is assumed 
that during phase change, these properties can be linearly inter- 
polated. It is proj >sed that Equation (2) is directly used to cal- 
culate the heat ol phase change. because this heat is evolved over 
a finite period of time, i.e., during the phase change, h rein'esents 
an ins tantaneous rate of heat dissipation. If this occurs over a 
time increment of AO ^ then from Etiiuition (2) 
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and the total heat oi pb.ase change is 
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Equation (48) is solved numerically to determine the heat of phase 
change of the sample media* 

5* SUMMARY 

It has been shown that important thermodynamic properties of a 
sample mass can be determined by comparison with samples of known 
thermodynamic properties such as distilled water* This is possible 
because of the unique capabilities of the device discussed in this 
report* A sensitivity analysis of the proposed equations remains to 
be accomplished. This may show that a more optimum approach sliould bo 
used. In any case* sufficient data are available for the determination 
of the thermodynamic properties of the sample mass. These properties 
may well hold an important key in the understanding of changes experi- 
enced by biological tissues during freeze, thaw or low temperature 
preservation. If such is the case, an understanding of these phenomena 
may lead to the development of optimum protocols for the frozen pre- 
servation of valuable biological tissues. 
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APPENDIX F 


RED BLOOD CELL TEST DATA 
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